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fH , In the Nambu-Jona-Lasinio model with Polyakov loops, we explore the relation between 

^ ^' the deconfinement and chiral phase transitions within the mean-field approximation. We 

I ' focus on the phase structure of the model and study the susceptibilities associated with 
corresponding order parameters. 

1. Introduction 

I Quantum Chromo dynamics (QCD) exhibits dynamical chiral symmetry breaking 

' and confinement. Both features are related with global symmetries of the QCD 

Lagrangian. While the chiral symmetry is exact in the limit of massless quarks 
the Z(3) center symmetry, which governs confinement, is exact in the opposite 
limit, i.e., for infinitely heavy quarks. The susceptibilities of the corresponding order 
parameters were calculated in lattice gauge theory (LGT) with finite quark masses^^ 
and hence both symmetries are explicitly broken. Nevertheless, the susceptibilities 
show a peak structure. Thus, the global symmetries of the Lagrangian are still 
relevant even for finite quark masses. Based on this consideration one can construct 
an effective Lagrangian in terms of the order parameters, the Ginzburg-Landau 
effective theory. 

Recently, gluon degrees of freedom were introduced in the Nambu-Jona-Lasinio 
(NJL) Lagrangian ^ throug h an effective gluon potential expressed in terms of 
Polyakov loops (PNJL model) SEl xhe model has a non- vanishing coupling of the 
constituent quarks to the Polyakov loop and mimics confinement in the sense that 
only three-quark states contribute to the thermodynamics in the low-temperature 
phase. Hence, the PNJL model yields a better description of QCD thermodynamics 
near the phase transition than the NJL model. Furthermore, due to the symmetries 
of the Lagrangian, the model belongs to the same universality class as that expected 
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for QCD. Thus, the PNJL model can be considered as a testing ground for the 
critical phenomena related to the breaking of the global ^(3) and chiral symmetries. 

It has been shown that the PNJL model, formulated at finite temperature and 
finite quark chemical potential, well reproduces some of the thermodynamical ob- 
servables calculated within LGT. The properties of the equation of state the 
in- medium modification of meson masses ^ as well as the validity and applicability 
of the Taylor expansion in quark chemica l potential used in LGT were recently 
addressed within the PNJL model I'^ I ^ I ^ I ^O I. In Ref. E] the model was extended to 
a system with finite isospin chemical potential and pion condensation was studied. 

Enhanced fluctuations are characteristic for phase transitions. Thus, the ex- 
ploration of fluctuations is a promising tool for probing the phase structure of a 
system. The phase boundaries can be identified by the response of the fluctua- 
tions to changes in the thermodynamic parameters. In this article we describe the 
susceptibilities of the order parameters and their properties in the PNJL model. 

2. Nambu Jona-Lasinio model v^rith Polyakov loops 

An extension of the NJL Lagrangian by coupling the quarks to a uniform temporal 
background gauge field, which manifests itself entirely in the Polyakov loop, has 
been propose d to account for interactions with the color gauge field in effective 
chiral models The PNJL Lagrangian is given by 

C = /:njl (V', V^, f [A] , $ [A] ; T, a*) - W(<i> [A] , $ [A] ; T) . (1) 

The interaction between the effective gluon field and the quarks in the PNJL La- 
grangian is implemented by means of a covariant derivative in £njl- The effective 
potential U of the gluon field in ([T|) is expressed in terms of the traced Polyakov 
loop $ and its conjugate $ 

$.i-Tr.L, $ = ^Tr.Xt, (2) 

where L is a matrix in color space related to the gauge field by 

r rfS 
L{x) — ■Pexp i I 
Jo 



(3) 



with V being the path (Euclidean time) ordering, and P = 1/T with A4 — zAq. In 
the heavy quark mass limit, QCD has the Z(3) center symmetry which is sponta- 
neously broken in the high-temperature phase. The thermal expectation value of 
the Polyakov loop ($) acts as an order parameter of the Z(3) symmetry. Conse- 
quently, (4>) = at low temperatures in the confined phase and ($} 7^ at high 
temperatures corresponding to the deconfined phase. 

The effective potential U{^, $) of the gluon field is expressed in terms of the 
Polyakov loops so as to preserve the Z{3) symmetry of the pure gauge theory tl^. 
We adopt an effective potential of the following form 
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Fig. 1. The chiral Xmm (dashed-line) and the Polyakov loop Xii (solid-line) susceptibilities in the 
chiral limit as functions of temperature T for fi = Q. 



with 

MT) = ao + ai(|) +«^(|)'+«3(|)'. (5) 

The coefficients Tq, Ui and bi are fixed by requiring that the equation of state 
obtained in pure gauge theory on the lattice is reproduced. In particular, at Tq — 270 
MeV the model reproduces the first order deconfinement phase transition of the 
pure gauge theory. 



3. Susceptibilities and the phase structure 

In the PNJL model the constituent quarks and the Polyakov loops are effective fields 
related with the order parameters for the chiral and Z{3) symmetry breaking. In 
LGT the susceptibilities associated with these fields show clear signals of the phase 
transitions ^. 

In Fig. [1] we show the chiral Xmm and Polyakov loop Xii — — (^) {^) suscep- 
tibilities computed at /x = in the PNJL model in the chiral limit within the mean 
field approximation. The chiral susceptibility exhibits a very narrow divergent peak 
at the chiral critical temperature Teh, while the peak of Xu is much broader and the 
susceptibility remains finite at all temperatures. This is due to the explicit breaking 
of the Z{3) symmetry by the presence of quark fields in the PNJL Lagrangian. Nev- 
ertheless, Xil still exhibits a peak structure that can be considered as a remnant for 
the deconfinement transition in this model. One finds interference of XiT and Xmm, 
as shown in Fig. [TJ At the chiral transition, T = Tch, there is a narrow structure 
in XiT- We stress that this feature is not necessarily related with deconfinement 
transition, but rather expresses the coupling of quarks to the Polyakov loops. Thus, 
for the parameters used in the model, the deconfinement transition, signaled by the 
broad peak in Xih sets in earlier than the chiral transition at vanishing net quark 
density. 

The peak positions of the Xmm and Xii susceptibilities can be used to determine 
the phase boundaries in the (T, /i)-plane. At finite chemical potential, there is a 
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Fig. 2. The phase diagram of the PNJL model in the chiral limit. The solid (dashed) line denotes 
the chiral (deconfinement) phase transition respectively. The TCP (bold-point) is located at (Tc = 
157, Ale = 266) MeV. 

shift of the chiral phase transition to lower temperatures. In Fig. [2] we show the 
resulting phase diagram for the PNJL model. The boundary lines of deconfinement 
and chiral symmetry restoration do not coincide. There is only one common point 
in the phase diagram where the two transitions appear simultaneously. Recent LGT 
results for two flavor both at vanishing and at finite quark chemical potential show 
that deconfinement and chiral symmetry restoration appear in QCD along the same 
critical line In general it is possible to choose the PNJL model parameters such 
that the both critical temperatures coincide at /i = 0. 

From Fig. [2] one finds that the slope of Tdcc as a function of /i is almost fiat, 
indicating that at low temperature the chiral phase transition should appear much 
earlier than deconfinement. However, there are general arguments, that the decon- 
finement transition should precede restoration of chiral symmetry (see e.g. 1131131) 
In view of this, it seems unlikely that at T ~ the chiral symmetry sets in at the 
lower baryon density than deconfinement. In the PNJL model, the parameters of 
the effective gluon potential were fixed by fitting quenched LGT calculations. Con- 
sequently, the parameters are taken as independent on /i. However, it is conceivable 
that the effect of dynamical quarks can modify the coefficients of this potential, 
thus resulting in /i~dependence of the parameters. The slope of Tdcc as a function 



of /i could be steeper Q Consequently, the effective Polyakov loop potential ([4]) with 
/i-indcpendent coefficients should be considered as a good approximation only for 



While the susceptibilities Xmm and Xii exhibit expected behaviors associated 
with the phase transitions, the diagonal Polyakov loop susceptibilities, 




A*/T < 1. 



(6) 



^ Such a modification was explored in Ref. 1151 where explicit fi- and Nf- dependence of To is 
extracted from the running coupling constant as , using arguments based on the renormalization 
group. 
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are negative in a broad temperature range above Tch This is in disagreement 
with recent lattice results, where in the presence of dynamical quarks xu is found 
to be always positive A possible cause for this behavior could be traced to the 
form of the effective Polyakov loop potential used in the Eq. (j4]). 

Recently an improved effective potential with temperature-dependent coeffi- 
cients has been suggested^ 



W($,l>;r) _ a{T) 



rp4 



+ b{T) In 



where 



a(T) = flo -I- ai 



(7) 



(8) 



The polynomial in $ and used in ([4]), is replaced by a logarithmic term, which 
accounts for the Haar measure in the group integral. The parameters in ([T]) were 
fixed to reproduce the lattice results for the thermodynamics of pure gauge theory 
and for the T-dependence of the Polyakov loop. The improved potential indeed 
yields positive values for all the Polyakov loop susceptibilities . We note that the 
phase diagram calculated with the improved potential is similar to that obtained 
with the previous choice of the Polyakov loop interactions, shown in Fig. [21 



4. Summary and discussions 

We introduced susceptibilities related with the three different order parameters in 
the PNJL model, and analyzed their properties and their behavior near the phase 
transitions. In particular, for the quark-antiquark and chiral density-density corre- 
lations we have discussed the interplay between the restoration of chiral symmetry 
and deconfinement. We observed that a coincidence of the deconfinement and chiral 
symmetry restoration is accidental in this model. 

We found that, within the mean field approximation and with the polynomial 
form of an effective gluon potential the correlations of the Polyakov loops in the 
quark-quark channel exhibit an unphysical behavior, being negative in a broad 
parameter range. This behavior was traced back to the parameterization of the 
Polyakov loop potential. We argue that the .Z(A^)-invariance of this potential and 
the fit to lattice thermodynamics in the pure gluon sector is not sufficient to provide 
a correct description of the Polyakov loop fluctuations. Actually it was pointed 
out ^ that the polynomial form used in this work does not possess the complete 
group structure of color SU(3) symmetry. The improved potential of Ref. 'Q'yields a 
positive, i.e. physical, xu susceptibility, in qualitative agreement with LGT results. 
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